Minimal model to describe the magnetism of CuGeOs 
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We show that to describe properly the low energy excitations of CuGeOa one must include the 
effects of a transverse antiferromagnetic coupling, which is estimated to be Jj_ ~ 0.15 J. Owing 
to this coupling the frustration in the chains is significantly lower than recent ID estimates based 
on purely one-dimensional arguments, we find J2 = 0.2J. Furthermore we have found a strong 
modulation of the nearest neighbour coupling due to the static distortion S — 0.065 , which is 5 
times higher than that previously deduced from a ID chain approach. Our set of parameters gives, 
i) a value of the distorsion which agrees well with some recent estimate for a lower bound, we are 
able to perfectly reproduce ii)the dispersions, iii) the experimental susceptibility at both high and 
low temperatures. By performing DMRG calculations for 2 coupled chains we have analysed the 
effect of the transverse coupling on the ratio of singlet to triplet gaps. The ratio is very sensitive 
to the parameters and the universality reported in the strict one dimensional case is lost. As an 
additional point, we provide a simple picture to explain the interesting new feature observed in 
recent inelastic neutron scattering experiments: the existence of a second branch of excitations. 



PACS numbers: 71.27.+a, 75.40. Mg, 75.90. +w 



I. INTRODUCTION 

The discovery of the non-organic compound CuGeOs 
has attracted considerable attention as a laboratory for 
low dimensional many-body quantum mechanics. This 
compound is believed to exhibit a spin Peierls transition 
at T = Tsp. Below this temperature the ground-state 
is dimerized, and simultaneously a gap opens in the ex- 
citation spectrum. Most of the attempts to extract the 
appropriate magnetic couplings have been in the frame- 
work of purely one dimensional system. This restriction 
to a one -dimensional picture was mainly justified by the 
features of Inelastic Neutron Scattering, the dispersion is 
the largest in the chain direction (c-direction) The 
model which has been widely studied that with dimer- 
ization and frustration; 

Hid ^JcY, ([1 + 5{-iy]S, ■ S,+i + aS, ■ S.+a) (1) 

i 

S measure the distorsion of the lattice at zero tem- 
perature and a is the measure of frustration, the ratio 
of second-nearest to nearest neighbour antiferromagnetic 
exchange. Attempts to fit the susceptibility data and 
the singlet-triplet gap have provided two different set of 
parameters {Jc,a,S): 

a) ( = 150K, a = 0.24, S = 0.030) in ref. (| 

b) (Jc = WOK, a = 0.36,(5 = 0.016) in ref. § The 
discrepancy in these two sets of parameters came from 
differing emphasis in the fitting procedures: In reference 
H the value a was first constrained to be at most 0.24 by 



the observation that the triplet gap scaled with temper- 
ature as a power of the lattice distortion and should thus 
be less than the (one-dimensional) values giving a spon- 
taneous dimerization. On the other hand it was made as 
large as possible to give a reasonable fit to the suscepti- 
bility. The distorsion S was fit from the dispersion at low 
temperatures. In reference Q the value of a was simply 
chosen to make the best fit possible to the susceptibility 
and S was taken to fit the triple gap at zero wave vector. 

Recently, Fabricius et al. [ij carefully reexamined the 
fit of the susceptibility to a strictly one-dimensional 
model. They deduced a value of a = 0.354, which is al- 
most identical to the estimate of Riera and Dobry. Their 
fit reproduced perfectly the susceptibility from 40^ up 
to lOOOX. For this reason the value a = 0.36 has been 
prefered to a = 0.24 and has been considered as the value 
from susceptibility measurements. 

In fact the issue is not settled, as apart from the ob- 
served relation of the gap to the lattice distortion, Ra- 
man scattering experiments lead to a value of a which 
is very much smaller, as we shall see: a value of a close 
to 0.2 was suggested. The argument for this is as fol- 
lows: It is known that the a system diescribed by (1) 
can exhibit a singlet bound-state excitation below the 
continuum A recent detailed analysis have shown 
that it appears at q — momentum when the frustra- 
tion is strictly non zero, for any given distorsion of the 
lattice 1^ . The most interesting feature is that the ratio 
R = As/ At of the gap to singlet excitations measured 
in Raman scattering and the gap to triplet excitations 
as visible in inelastic neutron scattering was shown to be 
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a universal function of the frustration parameter, i.e is 
independant of the dimerization parameter. This then, 
provides a direct tool to measure the amplitude of the 
frustration. Singlet excitations are directly observable in 
Raman scattering experiments. Such a bound state was 
first observed by Kuroe et al. 0]. The excitation energy 
is As ~ 1.78A( it could be directly concluded from 
fig.l in ref. §] that this leads to a « 0.20. Note that 
this value is very close to that proposed by Castilla et al 
In contrast a value of a « 0.35 would give R = 1.50, 
inconsistent with the experimental value. 

In conclusion, a description purely in terms of uncou- 
pled chains lead to a conflict: 

i) Susceptibility data ~* a — 0.36. 

ii) Raman scattering data ^ a « 0.2. 

The two values proposed suggest very different picture 
of the spin-Peierls transition: for a single chain there is 
a critical ac = 0.24 below which dimerization requires 
a coupling to phonons, and above which even without 
coupling to the lattice the magnetic chain spontaneously 
develops dimer order. While this may seem an academic 
point, given that a lattice dimerization is observed, but 
it is not if we wish to give a quantitative account of the 
physics: the estimate of that exchange dimerization de- 
pends on the frustration. It should be also stressed that 
Braden et al. have recently given an estimate of the lower 
bound value of the distorsion parameter wich has clearly 
excluded the value S = 0.016 of scenario b) Indeed 
they have found that S > 5min = 0.03. To solve this 
issue some have argued that an adiabatic approach was 
not appropriate and that the phonons degrees of freedom 
should be treated exactly ('dynamical phonon picture') 
. All recent attempts to include dynamical phonons 
were made in the framework of strictly ID system. Un- 
fortunately, because of the large size of the Hilbert space, 
the exact calculations are restricted to very small cluster 

In the following we will show that if we take into ac- 
count properly the transverse coupling between chains, 
these fundamental issues are solved. Within our ap- 
proach, described in the next section, we will be able 
to reproduce perfectly many experimental features and 
provide a more realistic value value of the distorsion pa- 
rameter 5. 



II. MIMIMAL 2D MODEL 

As was already mentioned in the first experimental 
papers, CuGeOs is not a strictly one-dimensional com- 
pound: a naive spin wave approach provide an estimate 
of the antiferromagnetic coupling in the perpendicular 
direction to the chain (b-direction) , Jf, w O.lJc, which is 
not so small 0Jl^. Unfortunately, the problem with a 
two-dimensional system is that much less is known both 
analytically and numerically. Numerical methods such as 
Exact Diagonalization and Density Matrix Renormaliza- 



tion Group can resolve only small clusters far from the 
thermodynamic limit. 

Most attempts to study CuGeOs as a quasi 2D com- 
pound were essentially within the framework of mean 
field theory |^]. An alternative, and powerful, approach 
for weakly coupled low-dimensional systems is to main- 
tain a precise treatment of the one dimension and treat 
only the interchain coupling perturbatively. We shall fol- 
low this method and supplement it with exact treatment 
of two chains for certain parameters as a check on the nu- 
merical parameters found. The low temperature phase of 
CuGeOs actually has a checker-board structure : the 
dimerization alternates from chain to chain. The mini- 
mal magnetic Hamiltonian is then: 

H2D = Hid + ^ ^i,r ■ Si,r+i (2) 

with, 

i,r 

(3) 

The variable r counts the different chains. The fac- 
tor (—1)'" explicitly takes into account the fact that the 
structure of the real system is checker-board. 

Starting from the limit of strong dimerization (in which 
triplet excitations are on nearest neighbour sites), we can 
treat the coupling between the chains perturbatively. It 
is easy to find, at the lowest order, that the dispersion is, 

uj{k\\,k±) =cji£i(fc||) - J^cos{k\\)cos{k_i) (4) 

where fcy is the wave- vector along the strongly coupled 
direction and k± is transverse. Eq. ^ implies the simple 
relation: 

A - Aid - J± (5) 

A is the real gap (including J±), and Aid denotes the 
gap in absence of coupling between the chains. 

Note that if the structure were not checker-board, but 
simply repeated in the b direction the dispersion would 
be different: u!{k\^,k±) = ujiD{k\\) — J±cos{k±). 

While equation (4) is strictly true for very strong 
dimerization, ie J± << w(fc||) it is a very good approxi- 
mation even when this is not satisfied, as we have verified 
by diagonalisation of two chains using Density Matrix 
Renormalization Group (DMRG). We then exploit this 
equation and use the exact one-dimensional dispersion 
from exact diagonalisation Ci;i£i(fc||). From exact diago- 
nalisation, in the parameter region relevant for our dis- 
cussion, a good fit for the dispersion is , 

^iD{k) — Jc\/ a — bcos{2k) (6) 

where 0=5 + ^ + and 6 = i + 
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III. FIXING THE PARAMETERS. 

In order to fix the parameters of the model, let us first 
start with the transverse coupling J^. For this purpose, 
we will use the experimental data for the dispersion as 
measured by inelastic neutron scattering . From equa- 
tion (^) the dispersion in the b-direction gives the trans- 
verse coupling from the difference beween the spin wave 
energy at qb the band edge and the centre, one has di- 
rectly 

i) = (5.6 - 2.0)/2 = l.^meV 

From equation (Q) the width of the ID dispersion, 
w(7r/2, A:^) — aj(0, k±) = a;iD(7r/2) — a;iD(0) = Aw. From 
eq. (H) the width of the ID dispersion is Aw « Jc. Thus 
using the given experimental data, we find, 

ii) ~ VlflmeV = \mK 

Hence, from i) and ii) we obtain Jj^ w 0. 15 Jc, this 
value of the transverse coupling confirms the poor ID 
character of CuGeOz- Note that the ratio of couplings is 
somewhat larger than estimated by the experimentalists 
who used an expression appropriate for a gap induced by 
anisotropy. As the measured anisotropy in spin is small 
p6t this is not appropriate. 

Now our task is to fix the value of the two remaining 
parameters, a and 5. As said before, this was easy in the 
ID picture, since the ratio singlet-triplet gap was shown 
to be a universal function of the frustration parameter 
only. Unfortunately, this is not the case anymore in the 
2D model {J± ^ 0). This will be illustrated later on. An 
alternative is to make use the high-temperature suscep- 
tibility data. As a reference for the high-T data, we have 
recalculated for the one-dimensional system the suscep- 
tibility with a — 0.35 and Jc = IbQK and g — 2.25, since 
with this set of parameters, the experimental data have 
been perfectly reproduced from -iOK up to lOOOK [Q. It 
is easy to get the first terms in the high T expansion for 
the susceptibility, in this limit one gets straigthforwardly, 

x{T)^^{9^iBf{a,-^) (7) 

where — j and ai — ^{\ + a + J^/ Jc) 

For the purely ID case, i.e J^ — 0, one requires a value 
of ai£) = 0.35. Thus the condition to reproduce the high 
T part of the susceptibility is, 

areel + "T^ = ~ 0-35 (8) 

Since Jj_ = 0.15Jc, the previous equation fixes unam- 
bigously the amplitude of the frustration: 

iii) ttrcel = 0.20 

From this relation, it can be concluded that a purely 
ID approach overestimates the real value of the frustra- 
tion parameter. In other words, to reproduce the High 
T behavior of the susceptibility in the one-dimensional 
limit, one requires a larger value of the frustration pa- 
rameter. 



As a check, we have performed the exact calculation of 
the susceptibility for a two chains system. This was done 
for a 2 X 8 system, using periodic boundary conditions 
in both directions, this means for the two chain problem 
changing J± to 2 J±. For this size we can fully diag- 
onalize the Hamiltonian in each subspace and calculate 
the thermodynamic functions. Since we are interested 
in the high temperature phase we have set (5 = 0. Note 
that, at high temperature the number of chains does not 
limit the accuracy of the calculation. In fig.^ we see that 
the agreement between the experimental data and the 
high-T expansion is very good for sufficiently large tem- 
peratures. We also see that the agreement with the ex- 
perimental data is very good down to 100 K. One should 
not pay too much attention to the fact that the pure one- 
dimensional approach was very good down to AbK. In 
our case, the finite size effect are definitely stronger: the 
chains are shorter and at low temperature the number of 
chains is important in the present case we have consid- 
ered only two chains. Nevertheless the agreement goes 
beyond the strict applicablity of the leading term. 

The last step consists in fixing the single remaining 
parameter, i.e the static distorsion amplitude. For this 
purpose, we had to perform a calculation of the gap in the 
one-dimensional case, with the frustration parameter set 
to a = 0.2. Using eq. ||, one immediately deduces Ai£) = 
0.322 Jc. Following the method of ref. |^ to extrapolate 
the data in the thermodynamic limit, we have obtained 
straightforwardly, 

iv) 6 w 0.065. 

We remark that this value is significantly larger that 
previously reported: it is almost 5 times larger than the 
value obtained in the purely one-dimensional approach 
1^ and closer to that estimated from the structure [||. 
At this point all the parameters of the model have been 
fixed from the experimental data. In order to test the 
validity of our set of parameters, we must confront it to 
different tests. This will be done in the following section. 

IV. PARAMETER TESTS. 

The realistic value of 5 we have found is already a first 
check of the validity of our set of parameters. We now 
proceed to a second test; let us check that we can accu- 
rately reproduce the dispersion data in both directions. 
Using eq. (||) and ^ we have performed the direct cal- 
culation of the dispersion in the chain and transverse 
directions. The plot of figj^ is rather convincing, the 
agreement between experimental and theoretical data is 
excellent. 

Now, as a third check of the validity of our approach, 
let us perform the calculation of the low temperature sus- 
ceptibility (T < Tsp) and compare it to available exper- 
imental data. For this purpose we assume that 5{T) = (5, 
which is reasonable except very close to the transition 
point . Following the method of ref. |jl^ the suscepti- 
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bility is accurately given by: 



whith, 



xiT) 



z{T) = 



z{T) 



l + 3z(T) 



(9) 



(10) 



The expression for the dispersion is given by equation 
(i. 

We have plotted in fig. M both the calculated suscep- 
tibility at low T and the experimental data . We ob- 
serve that the agreement is very good up to T w O.SSTsp. 
Note that between Q.'UhTsp and Tsp, one naturally ex- 
pects a deviation from the experimental data due to the 
sudden drop of 5{T) when approaching the transition 
temperature. As a remark one has to keep in mind that 
we did not use any fitting parameter. 

At this point, tt will be also an interesting point to esti- 
mate with our set of parameters the energy of the singlet 
bound state. As a preliminary step, let us first analyse 
the effect of the transverse coupling on the singlet-triplet 
ratio. Indeed, as said in the introduction part, the one 
dimensionnal calculation have shown that this ratio does 
only depend on the frustration only, leading in this case 
to a K, 0.2, to get the experimental value Rexp ~ 1-8. As 
we are now considering a weak, but non-zero transverse 
coupling we must know whether this ratio will change 
substantially with interchain exchange. To test this we 
have carried out Density Matrix Rcnormalization Group 
calculations on the system of two coupled chains. 



V. EFFECT OF Jx ON THE SINGLET-TRIPLET 
RATIO R. 

We have performed numerical calculations by apply- 
ing the DMRG method jl^ on the model defined by the 
Hamiltonian in Eq. (|^). Since this method is more ac- 
curate for systems with free ends, we will consider our 
two coupled chain model with open boundary conditions 
in the long direction. An unfavorable consequence of the 
DMRG is that the total momentum is not a good quan- 
tum number. As our aim is to determine the singlet and 
triplet energy gaps at fc = momentum we have used the 
spin reflection symmetry to rule out excited levels belong- 
ing to other k values . The singlet gap was calculated 
from the energy difference of the two lowest lying energy 
levels in the = spin sector with odd parity under 
spin reflection while the triplet gap was obtained from 
the lowest levels of the — Q and — 1 subspaces. 
Another unfavorable consequence of open boundary con- 
ditions, as in the nearest-neighbor valence-bound conflg- 
uration of the bilinear biquadratic model pH , is that free 
S ~ 1/2 spins remain at the ends of one of the chains, 
giving rise to a fourfold degenerate ground state. This is 



exactly the case for any flnite chain length for the extreme 
dimerization limit {5 — 1) and is true asymptotically for 
S < \. Since the extra degeneracies make the analysis 
of the spectrum more difficult we have removed the two 
outmost spins on one chain and set the end-coupling to 
J = 1 + 5 on the other chain. It worth mentioning that 
the total symmetry of all states becomes opposite under 
the spin reflection symmetry as it had been for the orig- 
inal problem. In most of the calculations we have used 
the more accurate version of DMRG, the so-called finite- 
lattice method to determine the energies more precisely. 
To further improve the efficiency of the calculations, we 
have also included the left-right reflection symmetry and 
all the investigated states were targeted independently. 
Since for finite dimerization the system is gapped, it was 
adequate to keep 100 — 200 block states to have the trun- 
cation error below 10^^ for chain up to 100 sites. The 
absolute error of our calculation was estimated by com- 
paring the ratio of the singlet and triplet energies (i?) 
obtained for the ladder model using the DMRG proce- 
dure at ckc = 0.24, Jj^ = to the exact value R = -/S 
given by the sine-gordon model. For stronger dimeriza- 
tion the deviation was 10~^ while for weaker dimerization 
it was found to be 10~^. 

The results for the ratio of the singlet and triplet en- 
ergy gaps (i?) as a function of Jx in the strong dimer- 
ization regime 5 — 0.2) and 0.2 < a < 0.26 values are 
plotted in fig. ||. It is seen from the figure that there ex- 
ists a finite region of J_l where the value of R is below 2, 
thus we have confirmed the the bound state is stable even 
in the 2D case for not too large interchain couplings. As a 
second result, we have found that the triplet gap is a lin- 
ear function of J±_ in agreement with eq. ^ In general the 
coefficient of this linear function should be different from 
unity and indeed this is apparent from the two chain re- 
sults. In order to estimate parameters, however, we take 
the linear coefficient as unity as the two chain estimate is 
too heavily biased by finite ( transverse ) size: doubling 
the transverse coupling to take into account the periodic 
boundary condition clearly oversetimates the correlations 
in the two transverse directions. The nonlinear behavior 
of R follows from the nonlinear dependence of the the 
singlet gap on Jj_. Taking the parameters for the frus- 
tration and R obtained from the ID chain calculations, 
namely a = 0.2, R = 1.80, one finds the same value of R 
for a different parameter set: a ~ 0.26 J_l ~ 0.12. 

In order to check the effect of dimerization in the case 
of 2D system we have performed the similar calculations, 
but for a weaker dimerization parameter 5 — 0.065. The 
result for i? as a function of Jj^ is also shown on the fig- 
ure. It is apparent from the curves that the R depends 
on (5, thus the ratio is no longer a universal function of 
the frustration parameter for agiven finite interchain cou- 
pling. On the other hand, we have found the triplet gap 
is linear function of the interchain coupling and the slope 
does not depend on the strength of the static distortion. 

It is clear that for weaker dimerization the ratio gets 
very sensitive to the interchain coupling. For a fixed frus- 
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tration, increasing interchain coupling leads to the dis- 
appearance of the bound state into the continuum. The 
fact that the ratio exceeds 2 means only that this state is 
no longer well defined. For instance, we find that when 
a = 0.2 the singlet bound-state disappear into the con- 
tinuum for Jperp > 0.1. While the value of the ratio of 
singlet to continuum can no longer be used to extract 
directily a value of the frustration, the very existence 
of the state gives a lower bound on the frustration. To 
conclude this section, one can say that it is extremely 
difficult from this two chain analyse to reproduce the ex- 
perimental value of the singlet-triplet ratio: one requires 
very precise values of the parameters of the model. 



VI. EXISTENCE OF A SECOND BRANCH IN INS 

Up to now we have assumed that the magnetic struc- 
ture is of equivalent but coupled chains at high temper- 
ature and below the spin-Peierls transition an identical 
dimerization along each chain in a checker-board struc- 
ture, ie the dimerization out of phase from chain to chain. 
There are two modifications to this view: one which we 
take to be unconfirmed, the other which we take to be 
reliable. 

The first is that X-ray analysis of the room tempera- 
ture structure have suggested a larger unit cell [ p2|j23| 
even above the spin-Peierls transition, in contrast to 
that found by neutrons Q . Re-examination by neutron 
diffraction |24| confirmed the earlier neutron results that 
these distortions are absent in the crystals studied by 
neutron diffraction, so while the differences in the ex- 
perimental structures remain to be fully elucidated, we 
shall assume that the structure at high temperatures is 
as taken. 

The second modification we take as necessary to take 
into account. Recent inelastic neutron results ||2^ on the 
dynamics have shown a new magnetic mode, that the 
experimentalists call an optic mode. We shall now show 
that the dispersion and amplitudes of this second mode 
are well explained by a slight refinement of the Hamil- 
tonian considered so far in supposing two inequivalcnt 
chains at low temperatures. As the numerical differences 
are small, this alters very little the dispersion and the 
susceptibility at low temperatures. We now give a sim- 
plified version based on the limit of strong dimerization. 

It is enough for our purpose to consider a system which 
consists of two inequivalent chains A and B coupled anti- 
ferromagnetically through . Wc first define the follow- 
ing set of parameters {J a = Jc + dJ, Sa = 5 + dS, a a = 
a + da) and {Jb = Jc — dJ, 6b = S — dS, ub = a — da) 
respectively for the chain A and B where it is assumed 
that dJ, dS and da are small. We also define 
and as the excitation energy for each chain, and 

^B{k\\) the associated eigenvectors . In the 
strong dimerization limit, it is relatively easy to per- 
form the calculation of the eignemodes. For a given ku, 



there are two possible modes and w_ corresponding 
to k± = 0, or TT. 



C^+(fc||) = C^(fc||) + ■^5Lu{k^l)^ + [Jj_COs(fc||)]2, (11) 



W_(fc||) = Uj{k\\) - y^(5w(fc||)2 + [Jj_COs(fc||)]2 (12) 

Their corresponding eigenstates are. 



Vl/+(fc||)=Vl/^(fc||) + (l- 



*_(fc||) = *A(fc||)-(l + 



(5ti;(fc|| 



J±cos{kH) 



Suj{kH 



J±cos{k\\) 



)^B{k\l), (13) 



)*B(fc||), (14) 



where uj = {uja + w_b)/2 and Sto — {lua — wb)/2. For 
simplicity we have assumed that we are only working in 
the vicinity of fc|| = or tt and that Sto <ti Jj_ cos (fc||), 
which is reasonable since the chains are only weakly in- 
equivalent. Now the crucial question is: Will we effec- 
tively observe both branches in INS experiment? 

To answer this question we have to calculate the ma- 
trix elements, F+ = (^'+ | S{k) \ *gs) and = | 
S{k) I 'i'as), where | ^'gs) is the Ground-state wave 
function and S{k) = Sf{k\\) + e*'=^''S'|(fc||). When per- 
forming the calculation at the lowest order one gets. 



foil 



Jk±b 



J±cos{k\\ 



dfoil 



Jk±b 



) (15) 



F- = /o(l - e''^' 7^^^) - rf/o(l + e^'^") (16) 

where /o = (*(fc||) | S'^(fc||) | *gs) for an isolated 
chain, considering the 'average' parameters (Jc, a), dfo 
corresponds to its variation. If we now fix for instance 
/c_L = then the spectral weight in both branches is. 



^+ = ll/o|P(l + 



2 J_LCOs(fc|| 



6uj{kH 



2Jj_cos(fc||) ' 



\dfo\ 



(5(jj(fc||) 
2J^cos(fc||) 



(17) 



(18) 



It is clear from these expressions that one will effec- 
tively observe two excitation branches, one with a sig- 
nificant spectral weight (the main branch) and another 
one with a relatively small one (secondary) , in agreement 
with the experimental observation. Typically the order 
of magnitude of the ratio of the spectral weight is. 



2 J_Lcos(fc|| ) 



(19) 



5 



This very simple calculation can be refined to make more 
quantitative comparison with the available experimental 
data. Let us repeat again that the previous simplified 
expressions are only valid in the vicinity of fcy = or 
TT and we are considering only two chains: ie k±_ is by 
definition or tt. 



VII. CONCLUSION 

In this paper we have proposed a minimal model to 
describe the magnetism of CuGeO^. It is shown un- 
ambigously that the compound is two-dimensional (the 
coupling in the third direction is very small) . We clearly 
identify the origin of the conflicts in the one-dimensional 
approach: i)small value of S and ii)large value of a es- 
timated from the susceptibilty data not consistent with 
the Raman experiments. The susceptibility previously 
calculated in a pure one-dimensional picture has in fact 
strongly overestimated the amplitude of the frustration. 
The smallness of S previously quoted, is simply due to un- 
derestimating the ID gap (see eq. ||) As a test of our set 
of parameters, we now provide a value of the distorsion 
which is more realistic, we are also able to reproduce per- 
fectly the high temperature (uniform phase) and low tem- 
perature (dimerized phase) behavior of the susceptibility 
without any additional fitting parameter. We could also 
perfectly reproduce the dispersion in both directions. As 
mentionned above, it is difficult to reproduce the singlet- 
triplet ratio from a study of two coupled chains, because 
of the extreme sensitivity of the ratio to the parame- 
ters. However, we believe that it would be interesting 
to perform the calculation of the Raman intensity, this 
will provide a measurable effect of the transverse cou- 
pling. We believe that this model with the parameters 
b) (Jc = U6K,a = 0.2,(5 = 0.065, Jf,/Jc = 0.15) should 
be the starting point for further studies on CuGeOs . For 
instance, the need to consider two inequivalent chains, as 
explained from observation of the extra branch in the in- 
elastic scattering, can be included by taking two slightly 
different values of the parameter S. 

As we have remarked, a value of a well below the crit- 
ical value for a chain implies that understanding of the 
dimerization inevitably involves the coupling to the lat- 
tice: the observed lattice distortion is not a secondary 
effect. The change of estimate implies that calculation 
of dilute phase diagram is quantitatively very different. 
While we consider the model with these parameters as 
the correct minimal model consistent with what is cur- 
rently known. 

Copper Germanatc is more correctly considered as a 
(spatially) anisotropic spin system rather than as a quasi- 
one-dimensional spin chain. The one-dimensional analy- 
ses that have been used until now therefore needed to be 
quantitatively modified. The apparent confiict between 
parameter sets previously obtained from different exper- 
iments reflects the simpliflcation. While we can give a 



satisfactory account to the experimental results to date, 
their remain aspects that require precision: 

• A more precise value of the interchain coupling re- 
quires a quantitative calculation in the anisotropic 
two dimensional spin system. Our simple estimate 
based on the strong dimer limit Jb/Jc = 0.15 may 
be modified by this. This also implies that a and S 
will be also affected. 

• The method proposed |6| to fix the frustration from 
the ratio of the singlet excitation seen in Raman 
to the triplet in neutron scattering is correct for a 
strictly one-dimensional model but must be modi- 
fied in the two dimensional case. We have verified 
the linearity of the triple gap predicted by the ex- 
treme dimer limit by means of a numerically exact 
two chain result. The existence and the position 
of a bound state in the Raman spectrum certainly 
constrain the parameters, but in a way which de- 
fies precise calculation at present. That is why, it 
may be useful to extend calculations of the Raman 
scattering to include the form of the continuum. 

• As yet there is no detailed understanding of 
anisotropic terms in the Hamiltonian. The bulk 
susceptibilities differ ||l8| but this has been inter- 
preted in terms of anisotropy in the tensor of g fac- 
tors. This may change with better understanding 
of the polarized neutron scattering cross-sections 
and/or more quantitative theories of the ESR [p6| . 

• The observation of the second inelastic branch in 
the neutron scattering indicates that the magnetic 
unit cell is doubled by a small inequivalence of ad- 
jacent chains in the low temperature phase. This is 
also seen in recent ESR experiments [^.It would 
be interesting to explore what physical parameter 
is inequivalent. 

VIII. ACKNOWLEDGEMENTS 

We are grateful to Jean-Paul Boucher, Emilio Lorenzo, 
Beatrice Grenier and Louis-Pierre Regnault for discus- 
sion of unpublished experimental data and for Olivier 
Cepas and Andrew Wills for discussions of the struc- 
tures and anisotropies. G. Bouzerar thanks Ph. Nozieres 
for hospitality at the Institut Laue Langevin. Ors Leg- 
eza was partially supported by the Hungarian Research 
Fund (OTKA) Grant No. 30173 and the French Ministry 
of Research and Technology Fund. 



6 



[1] M. Nishi, O. Fujita, and J. Akimitsu, Phys. Rev. B 50, 
6508 (1994); L.P. Regnault, M. Ain, B. Hennion, G. 
Dhalenne, and A. Revcolevschi, Phys. Rev. B 53, 5579 
(1996); 

[2] G. Castilla, S. Chakravarty, and V.J. Emery, Phys. Rev. 
Lett. 75, 1823 (1995). 

J. Riera and A. Dobry, Phys. Rev. B51, 16098 (1995). 
K. Fabricius, A. Klximper, U. Low, B. Biichner, T. 
Lorenz, G. Dhalenne, and A. Revcolevschi, Phys. Rev. 
B 57, 1102 (1998). 

B.S Shastry and B. Sutherland, Phys. Rev. Letter 47, 
964 (1981), G.S. Uhrig and H.J. Schulz, Phys. Rev. B 
54, R9624 (1996). 

G. Bouzerar, A. P. Kampf and G.L Japaridze, Phys. Rev. 
B 58, 3117 (1998). 

H. Kuroe et al. , Phys. Rev. B50, 16468 (1994). 
P.H.M. van Loosdrecht, J. P. Boucher, G. Martinez, G. 
Dhalenne and A. Revcolevschi, Phys. Rev. Lett. 76, 311 
(1996). 

M. Braden, G. Wilkendorf, J. Lorenzana, M. Ain, G. 
Mclntyre, M. Behruzi, G. Dhalenne, and A. Revcolevschi, 
Phys. Rev. B 57, 1105 (1996). 

B. Blichner, H. F eskhe, A. P. Kampf, and G. Wellein 
cond-mat 9806021 



[25] J. E. Lorenzo, L.P. Regnault, J. P. Boucher, B. Hennion, 
G. Dhalenne, and A. Revcolevschi, EuroPhys. Lett 45, 
619 (1999). 

[26] H. Nojiri, H. Ohta, S. Okubo O. Fujita, J. Akimitsu, M. 
Motokawa |cond-mat /9906074 , 



D. Augier, D. Poilblanc, E. Sorensen, L Affleck, Phys. 
Rev. B 58, 9110 (1998); G. Wellein, H. Fehske and A. P. 
Kampf Phys. Rev. Lett. 81, 3956 (1998). 
L.P. Regnault, M. Ain,B. Hennion, G. Dhalenne, and A. 
Revcolevschi, Physica B 213& 214, 278 (1995); M.C. 
Martin, G. Shirane, Y. Fujii, M. Nishi, O. Fujita, J. 
Akimitsu, M. Hase, and K. Uchinokura, Phys. Rev. B 
53, R14713 (1996). 

W. Brenig, Phys. Rev. B56, 14441 (1997); G.S. Uhrig, 
Phys. Rev. Lett. 79, 163 (1997); J.Zang, S. Chakravarty 
and A.R. Bishop, Phys. Rev. B55, R14705 (1997);R. 
Werner, C. Gros, Phys. Rev. B 57, 2897 (1998) 
K. Hirota, D.E. Cox, J.E. Lorenzo, G. Shirane, J.M. 
Tranqueda, M. Hase, K. Uchinokura, H. Kojima, Y. 
Shibuya, I. Tanaka, Phys. Rev. Lett. 73, 736 (1994). 
G. Bouzerar unpublished. 

T. M. BriU, J. P. Boucher, J. Voiron, G. Dhalenne, A. 
Revcolevschi, and J. P. Renard Phys. Rev. Lett. 73, 1545 
(1994). 

M. Troyer, H. Tsunetsugu and D. Wurtz , Phys. Rev. B 
50, 13515 (1994). 

M. Hase, I. Terasaki, K. Uchinokura, Phys. Rev. Lett. 70 
3651 (1993). 

5. R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. 
Rev. B 48, 10345 (1993). 

6. Legeza, J. Solyom Phys. Rev. B 56, 14449 (1997). 

L Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. 
Rev. Lett. 59, 799 (1987); Commun. Math. Phys. 115, 
477 (1988). 

M. Hidaka, M. Hatae, L Yamada,M. Nishi, J. Akimitsu, 

J. Phys. Condens. Matter. 9, 809 (1997). 

L Yamada, M. Hidaka, M. Nishi, J. Akimitsu, J. Mag. 

and Mag. Materialsl77-181, 703 (1998). 

M. Braden, E. Ressouche, B, Biichner, R. Kessler, G. 

Heger G. Dhalenne, and A. Revcolevschi, Phys. Rev. B 

57, 11497 (1998). 



FIG. 1. Dispersion in both directions. The symbols are 
experimental data from taken from ref. |l[. 



FIG. 2. Data for the effect of J± on the Singlet- Triplet ra- 
tio R. The distorsion 5 is fixed to 0.2 (a) and to 0.065 (b). We 
have considered several values of the frustration parameter a. 



FIG. 3. Susceptibility at High Temperature. Indi- 
rect Comparison between high-T expansion and experimen- 
tal data." Reference" was recalculated in the ID picture for a 
L = 16 sites system using the parameters of ref. M . 



FIG. 4. Susceptibility at Low Temperature. We compare 
the experimental data |l8[ with the analytical calculations. 
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